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^ . We present near-infrared X-band slit spectra of the nuclei of 25 Seyfert 2 
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CO index. A clear anti-correlation between the observed CO index and the 
nuclear K — L color is present, suggesting that a featureless hot dust continuum 



galaxies in the CfA and 12 /mi samples. The strength of the CO absorption 
features at 2.3-2.4 /mi produced by stars is measured in terms of a spectroscopic 



heated by an AGN contributes significantly to the observed i^-band fluxes in 
the nuclei of Seyfert 2 galaxies. After correction for this AGN contribution, we 
estimate nuclear stellar JT-band luminosities for all sources, and CO indices for 
sources with modestly large observed CO indices. The corrected CO indices for 
10 (=40%) Seyfert 2 nuclei are found to be as high as those observed in star- 
forming or elliptical (=spheroidal) galaxies. We combine the i^-band data with 
measurements of the L-band 3.3 /mi polycyclic aromatic hydrocarbon (PAH) 
emission feature, another powerful indicator for star-formation, and find that the 
3.3 /xm PAH to JT-band stellar luminosity ratios are substantially smaller than 
those of starburst galaxies. Our results suggest that the 3.3 /mi PAH emission 
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originates in the putative nuclear starbursts in the dusty tori surrounding the 
AGNs, because of its high surface brightness, whereas the i^-band CO absorption 
features detected at the nuclei are dominated by old bulge (=spheroid) stars, 
and thus may not be a powerful indicator for the nuclear starbursts. We see no 
clear difference in the strength of the CO absorption and PAH emission features 
between the CfA and 12 jum Seyfert 2s. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert - 
infrared: galaxies 

1. Introduction 

According to the current unification paradigm of Seyfert galaxies — the most numerous 
class of Active Galactic Nuclei (AGNs) in the local universe — , the two types of Seyfert 
galaxies, type 1 (which show broad optical emission lines) and type 2 (which do not), are 
both powered by a central mass-accreting supermassive black hole. The different emission 
properties are explained by the presence of a torus-shaped structure of molecular gas and 
dust surrounding the central supermassive black hole (the so-called "dusty torus") and the 
different viewing angles toward the torus (Antonucci 1993). Since the dusty torus is rich in 
molecular gas, it is a natural site for starbursts to occur (Fabian et al. 1998). Such a nuclear 
starburst in the torus could play an important role for AGN physics (Wada & Norman 2002), 
so that its observational understanding is very important. 

Attempts to put observational constraints on the nuclear starbursts in the torus in 
Seyfert 2 galaxies have been made extensively, from UV-optical (e.g., Gonzalez Delgado 
et al. 1998), through infrared (e.g., Oliva et al. 1999; Ivanov et al. 2000; Imanishi 2002; 
2003), to millimeter wavelength ranges (Kohno et al. 2002). Among these methods, infrared 
diagnostics have two major advantages. First, since nuclear starbursts in the torus are 
expected to have some amount of dust extinction, infrared observations are better suited to 
investigate their properties than those in the UV or optical, because of the reduced effects of 
dust extinction in the infrared. Second, compared to millimeter observations, a systematic 
study of the nuclear starbursts in a statistically meaningful number of Seyfert 2s is possible 
with a reasonable amount of telescope time. Imanishi (2003) obtained infrared L-band (2.8- 
4.1 /jm) slit spectra of 32 Seyfert 2s in the CfA (Huchra & Burg 1992) and 12 (Rush, 
Malkan, & Spinoglio 1993) samples. Utilizing the polycyclic aromatic hydrocarbon (PAH) 
emission feature at 3.3 fim as an excellent tool to disentangle star-formation from AGN 
activity (Imanishi 2002), the author discussed the properties of nuclear starbursts in Seyfert 
2s in a quantitative manner. 
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Nuclear starbursts in Seyfert 2s can also be investigated through infrared K-band (2.0- 
2.5 fj,m) spectroscopy. Since the CO absorption features at A rest = 2.3-2.4 /im in the rest- 
frame are produced by stars, and not by a pure AGN, these features can be a good tool to 
distinguish between stars and an AGN (Oliva et al. 1999; Ivanov et al. 2000). Among the 
32 Seyfert 2s studied by Imanishi (2003), 25 sources were observed with SpeX at the IRTF 
3m telescope (Rayner et al. 2003) (see Table 2 of Imanishi 2003). Since Spex enables to 
obtain K- and L-band spectra simultaneously, the stellar emission in Seyfert 2 nuclei can be 
investigated also from the iT-band CO absorption features, independently from the L-band 
PAH diagnostic. In this paper, we present Lf-band spectroscopic results and combine them 
with the L-band results to provide a better understanding of Seyfert 2 nuclei. H = 75 km 
s _1 Mpc -1 , f2 M = 0.3, and £7a = 0.7 are adopted throughout this paper. 



2. Targets, Observations and Data Analysis 

Detailed information and the observing log of the 25 Seyfert 2s are summarized in 
Table 1. The original sample selection criteria and observing details are found in Imanishi 
(2003). In short, Seyfert 2s in the CfA and 12 fim samples at z = 0.008-0.035 and at 
declination > —35° were selected. The CfA and 12 /im Seyfert 2 galaxies were selected 
based on their host galaxy magnitudes in the optical (Huchra & Burg 1992) and IRAS 12 
fim fluxes (Rush et al. 1993), respectively. The observations were made with SpeX at the 
IRTF 3m telescope, using the 1.9-4.2 /im cross-dispersed mode with a 1'.'6 wide slit (Rayner 
et al. 2003), which provides a spectral resolution of R ~ 500 in the Lf-band. The 1'.'6 aperture 
used corresponds to physical scales of between 250 pc (z = 0.008) and 1 kpc (z = 0.035), 
so our spectra are suitable to investigate the nuclear starbursts, rather than ring-shaped 
circumnuclear starbursts which are found typically at distances of ~1 kpc from the nuclei 
(Le Floc'h et al. 2001). The sky was clear and the seeing at K was / . / 6-0 / .'9 (full-width at 
half-maximum; FWHM) throughout the observations. 

The data reduction of the Lf-band spectra was made using IRAF 1 in the same way as 
the L-band data, as described by Imanishi (2003). The Seyfert 2 and standard star spectra 
were extracted by integrating signals over l'.'S-S'.'O along the slit, depending on the actual 
signal profiles. The Lf-band magnitudes of the standard stars in Table 1 were estimated from 
their F-band (0.6 fim) magnitudes, adopting V — K colors appropriate to the stellar types 
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of individual standard stars (Tokunaga 2000). The Seyfert 2 spectra were divided by those 
of the corresponding standard stars, and then multiplied by the spectra of blackbodies with 
temperatures corresponding to those of the individual standard stars (Table 1). Appropriate 
spectral binning was employed, particularly for faint sources. The intrinsic CO absorption 
of the adopted standard stars (A-, F-, G-type main sequence stars) is too weak (Arnaud, 
Gilmore, & Cameron 1989) to affect our main conclusions about the CO absorption strength 
of the Seyfert 2 nuclei. 



3. Results 

Figure 1 shows the flux-calibrated K-band slit spectra of the 25 Seyfert 2s in our sample. 
Some sources showed a larger scatter in the spectra over the A Q bs = 2.00-2.02 fj,m and 2.05- 
2.07 fj,m range in the observed frame than their neighboring wavelengths, due to the presence 
of the Earth's atmospheric absorption features. Since these data points may have bad effects 
on the continuum determination, with which to estimate the strength of the CO absorption 
features, we only used spectra at A Q bs = 2.07-2.5 fim. For some of the Seyfert 2 galaxies 
in our sample, spectra with a smaller wavelength coverage, in the close vicinity of the Brj 
emission line at A rcst = 2.166 ^m, are available (Veilleux, Goodrich, & Hill 1997; Rhee & 
Larkin 2000). For all the CfA Seyfert 2s in our sample except NGC 4501 and NGC 5347, 
Ivanov et al. (2000) presented moderate-resolution (R ~ 700) spectra at A rest = 2.0-2.4 fim 
and/or higher-resolution (R > 1000) spectra in narrower wavelength ranges near the CO 
absorption features. However, for the majority of the 25 Seyfert 2s, our spectra in Figure 1 
are the first moderate-resolution spectra with an almost complete iT-band coverage. Thus, a 
systematic and consistent investigation of the CO absorption features in a statistically larger 
number of Seyfert 2s is possible. 

To investigate the presence of nuclear starbursts in Seyfert 2s using the strength of 
the CO absorption features, we basically follow the method proposed by Doyon, Joseph, & 
Wright (1994), who defined a spectroscopic CO index as 

CO spcc = -2.51og< J R 2 . 36 >, (1) 

where <R 2 36 > is the average of actual signals at A rest = 2.31-2.40 /xm divided by a power-law 
continuum (F A = a x X 13 ) extrapolated from shorter wavelengths. When the CO absorption 
is present, <R2.36> should be less than unity and thus the CO spec should show a positive 
value. This CO index is insensitive to spectral resolution, because an average value within 
the wavelengths that show CO absorption is used. Thus, this index can be directly compared 
among independent spectra observed with slightly different resolutions. 
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To determine the continuum level, Doyon et al. (1994) used data points at A rcst = 2.0- 
2.29 jum, excluding clear emission lines, and extrapolated it to A rcst > 2.31 /im. In our 
sample, the redshifts are in the range z = 0.008-0.035 and only data points at wavelengths 
A bs — 2.07-2.5 /im are available. Thus, data points at wavelengths as short as A rcst = 2.0 
jum are not fully covered. We thus used data points at A rcst = 2.1-2.29 /im to determine the 
continuum level and extrapolated it to longer wavelengths. Following Doyon et al. (1994), 
a power-law shape was assumed for the continuum and clear emission lines, such as H 2 1-0 
S(l) (A rcst = 2.122 /im) or Br7 (A rcst = 2.166 /im), were excluded. Regarding the wavelength 
range that should contain the CO absorption features, our spectra of A Q b s = 2.07-2.5 /im 
covers A rcst = 2.31-2.4 /im for all the observed Seyfert 2s at z = 0.008-0.035, so that the same 
wavelength range of A rcs t = 2.31-2.4 /im as employed by Doyon et al. (1994) was adopted to 
measure the CO absorption strength. Although the wavelength range used for the continuum 
determination is slightly different, our spectroscopic CO index is essentially the same as that 
defined by Doyon et al. (1994). 

Whenever possible, we tried to use this method to determine the continuum level in 
a consistent way. However, in the case of NGC 262 and F04385— 0828, the extrapolated 
continuum determined in this way provided a negative value for CO spec , possibly because of 
the curvature of the continuum emission at A rcst = 2.1-2.29 /im. In these two sources, we 
used only data at A rcst = 2.22-2.29 /im and 2.18-2.29 /im, respectively, which provided more 
reasonable extrapolated continuum levels at A rcst = 2.31-2.40 /im. The adopted continuum 
levels are plotted as solid lines in Fig. 1. 

The values and histogram of the observed spectroscopic CO indices (C0 spcc _ o bs) are 
summarized in column 2 of Table 2 and Figure 2a, respectively. The largest uncertainty of 
the C0 spec _ o bs values comes from the ambiguity of the continuum determination, which is 
difficult to assess quantitatively. For some Seyfert 2s, we slightly changed the wavelength 
range used for the continuum determination as well as the removed data points which may 
be contaminated by emission lines. The C0 spcc _ o b s value differs by a factor of as high as 
0.03, which can be taken as a possible maximum uncertainty for C0 spcc _ o b s - 

For the majority of the CfA Seyfert 2s, Ivanov et al. (2000) estimated spectroscopic 
CO indices based on independent data. They used a very narrow wavelength range of 
A rest = 2.285±0.005 /im and 2.298±0.005 /im for the continuum and the data points which 
include CO absorption features, respectively. If this method were to be applied to our lower 
resolution spectra, the resulting uncertainty for C0 spcc _ o b s would be large, because only one 
or two data points would be included in each of these wavelength ranges. In Table 2 we 
compare the observed spectroscopic CO indices based on our wider wavelength coverage 
spectra, with the indices derived by Ivanov et al. (2000). For some sources (NGC 5929, 
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NGC 7674, and NGC 7682), both CO index measurements are in good agreement, but there 
is some discrepancy for other sources, possibly because of the different method to derive 
C0 spec _ o bs and/or different aperture size. 

In some objects, narrow (<500 km s _1 in FWHM) emission lines are discernible in the 
spectra. In Fig.l, we indicate the identification of some clear emission lines. In general, 
their equivalent widths are smaller than in starburst galaxies (Goldader et al. 1997b; Vanzi, 
Alonso-Herrero, & Rieke 1998; Coziol, Doyon, & Demers 2001). To quantitatively estimate 
these weak and narrow emission line fluxes, by tracing the line profiles, spectra with R > 
1000 are usually required. No quantitative estimates for the line fluxes are made based on 
our lower resolution (R ~ 500) spectra. 



4. Discussion 

4.1. CO Absorption Strength 

Unlike pure star-forming galaxies, a featureless continuum from hot dust heated by an 
AGN can contribute significantly to the observed X-band fluxes in Seyfert 2 galaxies. This 
contribution reduces the depth of the CO absorption feature; the degree of this reduction 
can vary, depending on the amount of AGN emission in the K-band. Indeed, the CO spe c-obs 
values of a large fraction of the observed Seyfert 2 galaxies in Fig. 2a are substantially smaller 
than the typical values observed in star-forming or elliptical (=spheroidal) galaxies (CO spec 
> 0.15; Goldader et al. 1997a,b; Ivanov et al. 2000). This is likely to be caused by dilution 
of the CO absorption feature by the AGN's featureless continuum emission. 

In the L-band, the emission from stars is much reduced when compared to that in the K- 
band, whereas hot (~1000K) dust heated by an AGN has strong emission both at K and L. 
Thus, the relative contribution from the AGN to the nuclear spectra increases substantially 
from K to L (Alonso-Herrero, Ward, & Kotilainen 1996; Alonso-Herrero et al. 2001). In 
fact, the typical K — L color of an AGN is usually larger (redder) than a normal stellar 
population (Willner et al. 1984; Alonso-Herrero et al. 2003). The K — L colors can be used 
to estimate the AGN contribution to the nuclear i^-band spectra. The nuclear K — L colors, 
derived from our slit spectra, are summarized in column 5 of Table 2. Figure 3a compares 
the nuclear K — L colors with the C0 spcc _ o b s values. The K — L color and CO spec -obs are 
anti-correlated in such a way that objects with redder K — L colors tend to show smaller 
C0 spcc _ o bs- The probability that a correlation is not present is found to be 0.01% (Isobe, 
Feigelson, & Nelson 1986). Since the observed K — L colors should become redder, when 
compared to those of a normal stellar population, for an increasing contribution from the 



-7- 



AGN, the anti-correlation in Fig. 3a suggests that the smaller C0 spec _ o bs values in objects 
with redder K — L colors are caused by the AGN contribution. This clearly demonstrates 
that to investigate the stellar emission from the CO absorption features in more detail, a 
correction for the AGN dilution is necessary. To make this correction quantitatively, we make 
use of the rest-frame equivalent width of the 3.3 PAH emission feature (EW3.3pah)- The 
EW3.3PAH values are summarized in column 6 of Table 2. 

Hot dust heated by an AGN produces a featureless L-band continuum, while star- 
forming galaxies usually show strong 3.3 ^m PAH emission with EW3.3PAH ~ 100 nm (Iman- 
ishi & Dudley 2000). Since the equivalent width is by definition robust to dust extinction, 
we assume that pure star-formation shows EW3.3PAH = 100 nm, while a pure AGN shows 
EW3.3PAH = nm. For the three Seyfert 2s with EW 3 . 3 pah values >40 nm (NGC 7682, Mrk 
938, and NGC 3660; Imanishi 2003), a significant contribution from star-formation to the L- 
band fluxes must be present. For these sources, we estimate the fractional stellar contribution 
at L, based on the observed EW3.3PAH values. More specifically, if the observed EW3.3PAH 
value is 70 nm, we regard that 70% of the observed L-band flux comes from star-formation, 
and the remaining 30% originates in AGN emission. We then extract L-band fluxes origi- 
nated in AGN activity from the observed L-band flux levels. For the Seyfert 2s with very 
small EW3.3PAH of < 20 nm (22 sources except NGC 7682, Mrk 938, and NGC 3660), we 
regard the observed L-band fluxes as dominated by AGN-powered emission (Imanishi 2003). 
Namely, the observed L-band fluxes are taken as AGN-originated L-band fluxes. 

To derive the AGN fraction in the Lf-band spectra from the estimated AGN-powered 
L-band flux levels, we assume that the intrinsic AGN colors in all Seyfert 2 galaxies are 
K — L = 2.79, as was done by Ivanov et al. (2000). Alonso-Herrero et al. (2003) have shown 
that the infrared AGN spectral energy distributions in Seyfert 2 galaxies are fitted with i u 
oc z/~ 3 , which provides K — L ~ 2.5-2.6. Thus, the adopted assumption of K — L = 2.79 
is reasonably justified from the observational point of view. After subtracting this AGN 
emission from the observed Lf-band spectra, we can infer the stellar Lf-band emission and 
thereby obtain corrected CO indices (CO spcc _ cor ). The derived stellar Lf-band luminosities 
are summarized in column 7 of Table 2. Regarding the corrected CO index (CO spcc _ cor ), if 
the AGN contribution to the observed Lf-band flux is estimated to be 60%, then the CO 
absorption optical depth increases by 2.5 {=1/(1—0.6)}, and the spectroscopic CO index 
increases accordingly. Thus, when the observed CO index is very small, say C0 spcc _ o bs 
< 0.1, due to dilution from AGN emission at K, a small uncertainty in the continuum 
determination can result in a large ambiguity in the corrected CO index. There are 10 
sources with C0 spcc _ o b s > 0.1. We regard the CO absorption as clearly detected in these 
10 sources, and estimate CO spcc -cor only for them, in order for the correction process to 
work reasonably. The CO spcc _ cor values of the 10 sources are summarized in column 3 of 
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Table 2, and the histogram of the CO indices corrected in this way is plotted in Figure 2b. 
Although the AGN correction certainly increases the CO indices, the change of their overall 
distribution is little. Thus, possible ambiguities of this AGN correction do not seriously 
affect our main conclusions. 

This method to correct the CO indices for AGN emission using our own slit spectra 
has important advantages over some other methods, such as that based on nuclear K- and 
L-band photometric data. First, since the K- and L-band spectra were taken simultaneously 
under the same observing conditions and probe the same regions in the observed Seyfert 2 
galaxies, there is little ambiguity on the flux inter-calibration. Second, the compact nuclear 
L-band emission observed in Seyferts is usually assumed to be dominated by AGN emission 
(Ivanov et al. 2000). Although this is true for the majority of Seyfert 2s (Imanishi 2002, 
2003), there are some exceptions. From our L-band slit spectra, we can clearly pick out these 
exceptions based on the EW3.3PAH values and take into account the stellar contamination to 
the nuclear L-band fluxes. 

However, there is one caution for this correction process. The PAH emission observed 
in star-forming regions primarily comes from photo-dissociation regions located between HII 
regions and molecular gas, and not from the HII regions themselves (Sellgren 1981). If 
the star-forming regions have such a high radiation density that their emitting volumes are 
predominantly occupied by HII regions with virtually no photo-dissociation regions, then no 
3.3 /iiii PAH emission is expected. Although this is the case for the nuclear star cluster in the 
dwarf galaxy NGC 5253 (Alonso-Herrero et al. 2004), this kind of star- forming regions are 
probably minor exceptions. Figure 3b compares the nuclear K — L color with EW3.3PAH for 
our sample of Seyfert 2s. Using a statistical test by Isobe et al. (1986), the probability that 
a correlation is not present is found to be 30%, suggesting an anti-correlation between the 
K—L and EW 3 3 pah, as expected from the assumption that the featureless L-band continuum 
originates in the AGN. If emission from pure HII regions were contributing significantly to 
the nuclear L-band featureless continuum in many Seyfert 2 galaxies, then the estimated 
AGN contribution to the nuclear Lf-band flux would be overestimated. In this case, the 
corrected CO indices would be smaller than the above estimates. 

In Fig. 2b, 40% (=10/25) of the observed Seyfert 2 galaxies have corrected CO indices 
(CO S p CC _ cor ) that are as large as those of star-forming or elliptical (=spheroidal) galaxies 
(CO S pcc > 0.15; Goldader et al. 1997a,b; Ivanov et al. 2000). Since very young star- 
formation, with less than a few million years in age, should show very small CO spcc values 
(Leitherer et al. 1999), such star-formation is missed in this Lf-band CO absorption method. 
Furthermore, when the AGN dilution is very severe in the Lf-band, the CO absorption 
features are very weak and so the correction of the observed spectroscopic CO indices for the 
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AGN dilution is not possible. In fact, for NGC 7674 and F01475-0740, although on-going 
nuclear star formation is indicated from the detected 3.3 //m PAH emission feature, no clear 
CO absorption features are detected. Thus, this fraction (=40%) should be taken as a lower 
limit of Seyfert 2 nuclei with stellar emission. 

4.2. The Nature of the Stellar Emission Detected in the K-band 

The stellar signatures detected in the nuclei of Seyfert 2s can be due either to nuclear 
starbursts in the AGN torus and/or old bulge (=spheroid) stars. Since both starbursts and 
old bulge stars show CO spcc > 0.15, these two populations are not easily distinguishable 
based on the CO spec values alone (Oliva et al. 1995). To further understand the nature of 
the detected stellar emission in the i^-band, we combine the K-band data with the L-band 
spectra (Imanishi 2003). Figure 4 compares observed CO index with EW3.3PAH (shown in 
the column 6 of Table 2). They are weakly correlated. 

Imanishi (2003) implicitly assumed that the detected 3.3 PAH emission originates 
in the putative nuclear starbursts in the AGN torus. Here we test this assumption in a 
quantitative manner based on the 3.3 /im PAH emission luminosity measured in our aper- 
tures, because starbursts and old bulge stars are distinguishable in terms of emission surface 
brightness. Heckman (2001) estimated a star-formation rate per unit area of 10" 1 M yr _1 
kpc~ 2 as the lower limit beyond which the superwind activity, a characteristic property of 
starburts, can occur. This limit corresponds to an infrared surface brightness of Sir ~ 2 
x 10 42 ergs s _1 kpc" 2 (Kennicutt 1998), or a 3.3 /im PAH emission surface brightness of 
S3.3PAH ~ 2 x 10 39 ergs s _1 kpc -2 (Imanishi 2002). Among the 25 Seyfert 2s studied in this 
paper, the 3.3 /im PAH emission was detected by Imanishi (2003) in the seven sources; NGC 
7674, NGC 7682, Mrk 938, F01475-0740, F04385-0828, NGC 3660, and NGC 4968. The 
surface brightnesses of the 3.3 /mi PAH emission are summarized in the last column of Table 
2. For all the PAH-detected sources, the observed surface brightnesses significantly exceed 
the threshold required for starbursts. Therefore, for the PAH-detected Seyfert 2 nuclei, we 
can conclude that the detected 3.3 fim PAH emission must come from a nuclear starburst, 
and cannot be of old stellar origin. 

Figure 5 compares the i^-band stellar luminosity with the 3.3 /im PAH emission lu- 
minosity from the Seyfert 2 nuclei measured from our slit spectra. In infrared-luminous 
starburst galaxies, the i^-band to infrared luminosity ratios are found to be Lk/Lir ~ 10~ L6 
(Goldader et al. 1997a). The 3.3 /im PAH to infrared luminosity ratios in starbursts are 
estimated to be L3.3PAH/L1R ~ 10~ 3 (Mouri et al. 1990; Satyapal et al. 1995; Imanishi 2002). 
Thus, the 3.3 /mi PAH to K-band luminosity ratios are expected to be L 3 3P ah/Lk ~ 10~ L4 , 
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if both luminosities trace starbursts. This value is shown as a solid line in Figure 5. The 
observed ratios in the Seyfert 2s in our sample are smaller, by a large factor, than the value 
expected for starbursts. If the starbursts are very young, with ages of less than a few million 
years, then Lk/Ljr decreases (Goldader et al. 1997a), and thus the discrepancy between the 
observed and expected L 3 . 3 pah/Lk ratios is even larger. 

Figure 5 suggests that in the nuclei of Seyfert 2s, compared to starbursts, L 3 3P ah is 
underluminous or Lk is overluminous. One possibility for the L 3 . 3 pah underluminosity is the 
destruction of PAH emitting material by energetic radiation from the central AGN (Voit 
1992). However, since the putative nuclear starbursts are expected to occur at the outer 
part of the obscuring torus where the AGN radiation is significantly attenuated (Heckman 
et al. 1997; Imanishi 2003), this destruction is unlikely to be severe. Aside from the PAH 
destruction by the AGN radiation, the PAH emission may be suppressed if the stellar radia- 
tion density is so high that the volume fraction of photo-dissociation regions relative to HII 
regions is substantially smaller than typical starbursts. If the PAH destruction in the nu- 
clear starbursts in the AGN tori were more severe than typical starburst galaxies, the actual 
infrared luminosity of the nuclear starbursts would be substantially larger than the estimate 
based on the assumption of L 33 pah/Lir ~ 10~ 3 (Imanishi 2003). However, we regard this 
possibility as unlikely, because for a few selected Seyfert 2s for which the luminosities of the 
nuclear starbursts were estimated both from the UV and 3.3 fim PAH emission (assuming 
L 3 . 3 pah/Lir ~ 10~ 3 ), the estimated nuclear starburst luminosities agree quantitatively well 
(Imanishi 2002). 

We suggest that our K- and L-band results are more naturally explained by a sce- 
nario where the 3.3 ^m PAH emission originates in the spatially-unresolved (smaller than 
sub-arcsec) nuclear starbursts in the AGN torus, whereas the signatures of stellar emission 
detected in the i^-band spectra are produced mainly by old bulge stars. In fact, the nuclear 
L-band emission of Seyfert 2s is dominated by spatially-unresolved emission rather than the 
bulge component, whereas the extended stellar components often contribute predominantly 
to the nuclear K-band emission (Alonso-Herrero, Ward, & Kotilainen 1996; Alonso-Herrero 
et al. 1998, 2001). Old bulge stars can contribute significantly to the nuclear X-band fluxes, 
but cannot to the PAH emission due to the paucity of PAH-exciting far-UV photons in bulges 
compared to starbursts. This scenario explains the small L PA h/Lk ratios observed in Figure 
5. In this scenario, the CO absorption strength is primarily due to old bulge stars, rather 
than nuclear starbursts, and thus it may not be a useful tool to detect the putative nuclear 
starbursts in the AGN torus. This scenario is consistent with the conclusion by Ivanov et 
al. (2000). 
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4.3. Do the CfA and 12 /im Seyfert 2 Galaxies Contain Different Nuclear 

Stellar Properties? 

To improve the statistics, so far we have discussed near-infrared emission properties by 
combining both the CfA and 12 //m Seyfert 2s together. However, since selection criteria of 
these samples are different, a separate discussion of Seyfert 2s in each sample is of interest. 

Figure 6 plots the distribution of CO spcc _ obs , EW 3 . 3 pah, L 3 . 3 pah/Lir, and L 3 . 3 p A H/L K stciiar- 
In all the plots, no systematic differences are found between Seyfert 2s in the CfA and 12 jum 
samples. Thus, the main arguments in this paper are taken to be applicable to Seyfert 2 stel- 
lar populations in general. In particular, the L 3 3 pah/Lir ratios in Fig. 6c and L 3 3 pAH/LKsteiiar 
in Fig.6d trace the relative contribution from the nuclear starbursts to the infrared and nu- 
clear stellar K-band luminosity, respectively (Imanishi 2003). There was a possibility that 
the 12 //m Seyfert 2s may be biased to those with powerful nuclear starbursts (Ho & Ulvestad 
2001), because the starbursts can be strong emitters at 12yum. The absence of systematic 
differences in Fig.6c and 6d suggests that no such strong bias is present in the 12 /xm Seyfert 
2 sample. 

5. Conclusion 

We have presented results from a near-infrared iCband slit spectroscopic survey of 25 
Seyfert 2 galaxies in the CfA and 12 /im samples. These data have been combined with results 
from the 3.3 /im PAH emission survey of Imanishi (2003). The CO absorption features at 
A res t = 2.3-2.4 fim have been used to investigate the properties of the stellar populations in 
the nuclei of Seyfert 2 galaxies. We have found that the observed values of the CO indices 
in a large fraction of Seyfert 2 nuclei are affected by dilution produced by the featureless 
continuum associated with the AGN, because Seyfert 2 nuclei with redder K — L colors 
tend to show shallower CO absorption features. For sources with modestly large observed 
CO indices (>0.1), we have corrected the CO indices for AGN dilution. This correction has 
made a small change for the overall distribution of the CO indices in our sample and we 
have estimated that ~40% (=10/25) of Seyfert 2 nuclei have CO indices similar to those 
of star-forming or elliptical (=spheroidal) galaxies. This is only a lower limit of Seyfert 2 
nuclei with stellar emission because young star-formation shows no CO absorption features 
and because in Seyfert 2 nuclei where the iT-band light is dominated by AGN emission, the 
correction of the observed CO indices for AGN dilution is not possible. 

A comparison between the nuclear 3.3 /im PAH to stellar iCband luminosity ratio in 
the nuclei of Seyfert 2s and observations for infrared luminous starburst galaxies implies that 
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the CO absorption features detected in the nuclear K-band spectra of Seyfert 2 galaxies are 
due mainly to old bulge (=spheroid) stars rather than the putative nuclear starbursts in 
the AGN torus, and thus that the putative starbursts are difficult to detect with i^-band 
spectra over the scales probed by our data (250 pc to lkpc). We found no clear difference 
in the properties of the 2.3 /xm CO absorption and 3.3 /j,m PAH emission features between 
the CfA and 12 jum Seyfert 2 samples. 
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Note. — Column (1): Object. Column (2): Redshift. Column (3): Physical scale of the employed slit width of 1'.'6. Column (4): Observing date in 
UT. Column (5): Net on-source integration time in min. Column (6): Standard star name. Column (7): Adopted i\-band magnitude for the standard 
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NGC 4501 is classified as a LINER in the CfA sample (Huchra & Burg 1992). 
Flux calibration was made using HR 8457. 



Table 2. Near-infrared Properties of the Observed Seyfert 2 Galaxies 
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Fig. 1. — Infrared 2.07-2.5 /mi spectra of the 25 Seyfert 2 nuclei. The abscissa and ordinate 
are the observed wavelength in /im and flux F\ in 10~ 15 W m~ 2 /im" 1 , respectively. The solid 
lines are the adopted continuum levels, with respect to which the CO indices are measured. 
The arrows indicate the CO absorption features. When an observed CO index is > 0.1 (< 
0.1), the term "CO absorption" ["CO absorption (?)"] is used. Some strong emission lines 
are indicated; H2 1-0 S(l) emission line at A rest = 2.122 /im (denoted as S(l)), Br7 A rcst = 
2.166 /mi (Br 7 ), H 2 1-0 S(0) at A rcst = 2.223 /mi (S(0)), H 2 1-0 Q(l) at A rcst = 2.407 /im 
(Q(l)), and H 2 1-0 Q(3) at A rcst = 2.424 /im (Q(3)). 
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Fig. 2. — (^aj: Histogram of the observed CO index (C0 spcc _ o b s )- (b): Histogram of the 
corrected CO index (CO spcc _cor)- For sources with observed CO indices < 0.10, no AGN 
correction was attempted. Star-forming or elliptical galaxies typically show CO spec > 0.15 
(Goldader et al. 1997a,b; Ivanov et al. 2000). 
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Fig. 3. — (a): Nuclear K — L magnitude (abscissa) and observed spectroscopic CO index 
(ordinate). The uncertainty of the K — L color is <0.1 mag. (b): The ordinate is the rest- 
frame equivalent width of the 3.3 ^m PAH emission feature. Its uncertainty is at a level of 
<30% for the PAH-detected sources. 
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Fig. 4. — Observed spectroscopic CO index (abscissa) and rest-frame equivalent width of 
the 3.3 PAH emission feature (ordinate). 
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Fig. 5. — The observed nuclear stellar if-band luminosity, defined as AL A at 2.2 fim (abscissa) 
and nuclear 3.3 /zm PAH luminosity estimated from our slit spectra (ordinate). The errors 
of the PAH luminosities for the PAH-detected sources are at a level of <30%. The solid line 
indicates the correlation for infrared-luminous starburst galaxies (see text). 



-28- 



d 



d 



6 Csl 

9- o 



1 1 

(a) CfA 


■ 1 i 

CfA 
12/im 


i i i i | i i 

12yLtm 


o 
o 


■ i 
(b) CfA 


1 1 i 

CfA 
^2/J,m 


. . . . i . . 

Mfim 

* 






★ 


T 

c 


• 




★ 


• 

- • 


■ . 
■ ■ 


★ * * - 

★ 

* * *: 

★ * ★ 


X 

< 

to ^ 
to '~~ 

UJ 


: j 


■ 

I 


' * t 




. . T 








i 








(c) CfA 


1 1 I 

CfA 


,,,,,,, 


to 
i 




(d) CfA 


■ 1 i 

CfA 


i i i i i i i 

12//m 






12/zm 




o 


O 
CM 




12/zm 


★ 




■ 't 

• 


T 


★ : 


X 
n 


O 




I 


** * ** 1 


d 


■ 


^^'^^ '^^ r ^^'^^ r 

Y "v" v* ▼ "v* 

★ 


L K stell. 


ID 


: I 

:»* * 








T 


* * * ; 


\ 

X 

< 

CL 

to 
to 


CM 


** 












_J 










o 




i 


■ ■ ■ ■ i ■ ■ 




LO 




i 





Fig. 6. — (a): The distribution of the observed spectroscopic CO index. Solid circles: CfA 
Seyfert 2s. Solid stars: 12 Seyfert 2s. Five sources listed in both samples are plotted 
as solid squares, (b): The rest-frame equivalent width of the 3.3 /im PAH emission feature. 
(c): The 3.3 /im PAH to infrared luminosity ratio, taken from Imanishi (2003). (d): The 3.3 
^m PAH to stellar if-band luminosity ratio. 



